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CURVATURE-BASED FEATURE SELECTION FOR 
ELECTROPHYSIOLOGIC SIGNALS 

The present invention relates generally to cardiac rhythm management 
systems, and more particularly, it pertains to a system and method of 
classification and detection of cardiac signals. 

Background of the Invention 

Analysis of cardiac signals, which is routinely performed in 
electrocardiography is generally based on visual inspection to quantify or qualify 
wave morphology for the purpose of identifying and classifying abnormal 
patterns. Certain morphological characteristics of commonly recorded signals 
have high diagnostic value. The shape and inter arrival times of R-waves 
recorded in the electrocardiogram generally provide a wealth of information 
about the state of the heart. Accordingly, automated approaches for identifying 
and classifying abnormalities in signals such as cardiac signals have sought to 
use a signal's significant morphologic characteristics. 

However, given the wide diversity of possible shapes for cardiac signals, 
it is usually not possible for an automatic approach to identify significant 
characteristics that can be used for unambiguous classification. Rather, 
automated classification approaches generally compare the entire morphological 
shape of a signal with the shape of similar signals with known abnormalities but 
without particular regard to the specific characteristics that the signals contain. 
Alternatively, automated classification approaches restrict the automated 
examination only to those signals which are essentially normal and use detailed 
metrics (for example QRS width, QT interval or ST segment amplitude) of the 
essentially normal morphology for classifying abnormalities. 

Despite its importance in the analysis of biologic signals, the automated 
and accurate identification and quantification of the significant morphological 
characteristics (for example turns, peaks, knees, inflection points, and the like) in 
any cardiac signal (both abnormal as well as normal) is still in a developing 
stage. Existing methods have used the concept of sharpness (for example to 
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detect R-waves) but have had limited success. This is due in part to the overly 
simplistic mathematical treatment this concept has received, as reflected in the 
rudimentary algorithms used for these measurements. Most of the current 
detection methods rely on three point interpolations to measure sharpness. The 
simplest and most commonly used methods for measuring peaks of R-waves are 
based upon Taylor-series approximations to estimate the second derivative of the 
sensed signal. This formula utilizes highly local information (the point at the 
peak and its two close neighbors) ignoring nearby points which may contribute 
to signal peak. Other popular approaches utilize less local data, such as the peak 
and two adjacent extrema. All of these methods, which rely on three-point 
estimates of sharpness, may produce inaccurate estimates, if waveforms are 
complex or are contaminated with noise. Thus, a need exists for automated 
identification and classification of peaks, knees, inflection points, and the like in 
sensed cardiac signals that takes into account wave scale and complexity that can 
yield a more accurate estimate of peaks for identifying and classifying 
abnormalities in cardiac signals. 

Summary of the Invention 
The present subject matter provides a curvature based method of 
selecting features from electrophysiologic signals for purpose of complex 
identification and classification. According to one aspect of the present subject 
matter, this is accomplished by sensing a cardiac signal (sensing the cardiac 
signal includes sensing complexes continuously on a real-time basis) and 

computing curvatures on a sample point-by-sample point basis (X ls X 2 , X 3i 

Xj ) on the sensed cardiac signal on a continuous basis. In one embodiment, the 

curvature at the sample points X lf X 2 , X 3 X r are computed by fitting a 

cubic least square error curve (using N number of sample points) to the sensed 
cardiac signal. In this embodiment, N is an odd number, and the sample point 
(where the curvature is computed) is at a mid point of the N number of sample 
points. 

Also on a continuous basis, features of significant interest are extracted 
from the computed curvatures. In some embodiments, this is accomplished by 
comparing the computed curvatures at the sample points X u X2, X 3t X t toa 
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set of predetermined threshold values. In some embodiments, the features are 
extracted based on computing features such as a time when the feature occurs, an 
amplitude of the feature, and other similar features at each of the sample points 
and comparing them to a set of threshold values. 

Also on a continuous basis a set of features associated with a first 
complex in the sensed cardiac signal are identified and separated from the 
continuously computed and extracted features upon detecting a second 
subsequent complex. The second subsequent complex is a complex that is 
adjacent to the first complex and occurs substantially immediately after the first 
complex. This process of identifying and separating extracted features repeats 
itself from one sensed complex to another subsequent sensed complex on a real 
time basis. One reason for identifying and separating the set of features 
associated with the first complex is to prevent the features associated with the 
first complex from mixing with the features associated the second subsequent 
complex. Separating the features associated with first complex aids in classifying 
the sensed first complex. 

Next, the process includes identifying a fiducial feature from the 
separated set of extracted features associated with the first complex and aligning 
the separated set of features with respect to the identified fiducial feature. In one 
embodiment, fiducial feature is identified based on comparing the times when 
each of the separated set of features occur with a time when a complex 
associated with the separated set of features is detected on the sensed cardiac 
signal, and selecting a feature from the separated set of features that is closest in 
time to flie time when the complex associated with the separated set of features 
' was detected. One reason for using a time when a complex is detected (such as R 
wave) in identifying the fiducial feature, is because the detection of a complex in 
a sensed cardiac signal is generally more reliable and consistent. 

Next, the process includes aligning the separated set of features around 
the identified fiducial feature. Aligning the separated set of features around the 
identified fiducial feature aids in normalizing each of the separated set of feature 
around a datum such as the associated detected complex, and further aids in 
comparing the separated set of features with a set of predetermined templates. 



3 



WO 02/40094 PCT/US01/45763 

Next the process includes comparing the aligned set of features to a set of 
predetermined templates to classify the associated complex. In some 
embodiments, the predetermined templates are a set of identified complexes 
associated with known cardiac arrhythmias that would assist in comparing and 
classifying the extracted set of features. In some embodiments, a therapy is 
provided to the heart based on the outcome of the classification. The above 
described process repeats itself on a continuous basis for a real-time 
classification of complexes from the sensed cardiac signal. 

Other aspects of the invention will be apparent on reading the following 
detailed description of the invention and viewing the drawings that form a part 
thereof. 

RHftf Descri ption of the Drawings 

Figure 1 is a cardiac signal according to one aspect of the present invention. 

Figure 2 is a curvature figure derived from the cardiac signal of Figure 1 
according to one aspect of the present invention. 

Figures 3A and 3B illustrate separating and aligning features derived 
from the curvature figure of Figure 2 according to one aspect of the present 
invention. 

Figure 4 is a flow diagram illustrating generally one embodiment of 
operation of the present subject matter. 

Figure 5 is a schematic/block diagram illustrating generally one 
embodiment of portions of a cardiac rhythm management system of the present 
subject matter. 

Detailed Description 

In the following detailed description, reference is made to the 
accompanying drawings which form a part hereof; and in which is shown by way 
of illustration specific embodiments in which the invention may be practiced. 
These embodiments are described in sufficient detail to enable those skilled in 
the art to practice the invention, and it is to be understood that the embodiments 
may be combined, or that other embodiments may be utilized and that structural, 
logical and electrical changes may be made without departing from the spirit and 
scope of the present invention. The following detailed description is, therefore, 
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not to be taken in a limiting sense, and the scope of the present invention is 
defined by the appended claims and their equivalents. 

Analysis of cardiac signals, which is routinely performed in 
electrocardiography is generally based on visual inspection to quantify or qualify 
wave morphology for the purpose of identifying and classifying abnormal 
patterns. Certain morphological characteristics of commonly recorded signals 
have high diagnostic value. The shape and inter arrival times of R-waves 
recorded in the electrocardiogram generally provide a wealth of information 
about the state of the heart. Accordingly, automated approaches for identifying 
and classifying abnormalities in signals such as cardiac signals have sought to 
use a signal's significant morphologic characteristics. 

However, given the wide diversity of possible shapes for cardiac signals, 
it is usually not possible for an automatic approach to identify significant 
characteristics that can be used for unambiguous classification. Rather, 
automated classification approaches generally compare the entire morphological 
shape of a signal with the shape of similar signals with known abnormalities but 
without particular regard to the specific characteristics that the signals contain. 
Alternatively, automated classification approaches restrict the automated 
examination only to those signals which are essentially normal and use detailed 
metrics (for example QRS width, QT interval or ST segment amplitude) of the 
essentially normal morphology for classifying abnormalities. 

Despite its importance in the analysis of biologic signals, the automated 
and accurate identification and quantification of the significant morphological 
characteristics (for example turns, peaks, knees, inflection points, and the like) in 
any cardiac signal (both abnormal as well as normal) is still in a developing 
stage. Existing methods have used the concept of sharpness (for example to 
detect R-waves) but have had limited success. This is due in part to the overly 
simplistic mathematical treatment this concept has received, as reflected in the 
rudimentary algorithms used for these measurements. Most of the current 
detection methods rely on three point interpolations to measure sharpness, may 
produce inaccurate estimates, if waveforms are complex or are contaminated 
with noise. 
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General S ystem Overview 
The present subject matter provides, among other things, a system of 
classifying a cardiac signal by computing a curvature versus time signal from a 
sensed cardiac signal and comparing the computed curvature versus time signal 
with a set of predetermined templates. The present system also provides an 
automated identification and classification of peaks, knees, inflection points, and 
the like in cardiac signals that takes into account wave scale and complexity that 
can yield a more accurate classification of cardiac signals. 

The first step in the process includes computing continuously a curvature 
versus time signal from a sensed cardiac signal on a real-time basis. The process 
also includes extracting features (features are generally significant points of 
interest on the sensed cardiac signal) continuously on a real-time basis from the 
computed curvature versus time signal. Next, the process includes continuously 
separating a set of features associated with a first complex upon detecting a 
second subsequent complex so that each individual complex can be identified 
and classified on a real-time basis. Then the process includes aligning and 
classifying the extracted and separated set of features by comparing the separated 
set of features with a set of predetermined templates. Other aspects of the 
invention will be apparent on reading the following detailed description of the 
invention and viewing the drawings that form a part thereof. 

Referring now to Figure 1, there is one embodiment of the first step in 
the process of sensing a cardiac signal 100 according to the present invention. 
Shown in Figure 1 are detection of heart beats associated with first and second 
complexes 1 10 and 120 at times '// and %' associated with the sensed cardiac 
. signal 100 on a real-time basis on a time line 130. Also, shown on the time line 

130 are the various sample points X„ X 2 , X 3 X, used in computing a 

curvature versus time signal on the first complex 1 10. In the example 
embodiment shown in Figure 1, the sample points can be some fixed interval of 
time, such as 2, 4, 6... milliseconds, at which curvatures are computed on the 
sensed cardiac signal. Also shown in this embodiment, is the detection of the 
complexes (1 10, 120, and so on) on a continuous basis on the sensed cardiac 
signal 100. 
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Referring now to Figure 2, there is shown one embodiment of generating 
a curvature versus time signal 200 on a real-time basis derived from the sensed 
cardiac signal 100 shown in Figure 1 according to the present invention. Figure 2 
is also on the same time line 130 as in Figure L Curvature figure 200 shown in 
Figure 2 is a computed curvature (expressed in time) versus time signal 
generated continuously by computing curvatures at the sample points X l9 X 2 , 
X 3 Xj on the sensed cardiac signal 100 shown in Figure 1 . Curvature figures 
210 and 220 are computed curvature versus time signals associated with first and 
second complexes 1 10 and 120 shown in Figure 1, respectively. In this example 
embodiment, the curvature versus time signal is computed by overcoming the 
dimensionality of the curvature by changing the dimensionality of the curvature 
to time. The reason for overcoming the dimensionality of the curvature is to 
convert the time versus voltage signal, which is time versus curvature, to a time 
versus time signal to facilitate easier computation of curvature (for example, 
when a curve exists in length-length space, the quantity has dimension of length 
and is traditionally called the radius of curvature, and if the curve exists in time- 
time space, the analogous quantity will be called the radius of time). In this time- 
time space the curvature has a dimension of 1/Time. Also, shown in Figure 2 are 
some example sample points 230 at which curvatures are computed. 

Referring now to Figures 3A and 3B, there is shown one example 
embodiment of extracting and aligning 300 sets of features associated with 
respective detected complexes. Figure 3A, illustrates continuously extracting 
features 320 from the computed curvature versus time signal 200 shown in 
Figure 2. Also shown in this example embodiment, are some sample times in 
milliseconds 330 when the extracted features occur on the time line 130. Figure 
3A is also on the same time as that of Figures 1 and 2. In the example 
embodiment shown in Figure 3 A, the sample times (where the curvatures are 
computed) have a fixed interval of 2 milliseconds. Also, in this example 
embodiment, the beat associated with first complex 1 10 is detected at 1 1 
milliseconds (*,) and the beat associated with the second subsequent complex 
120 is detected at 101 milliseconds (tj. 
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Figure 3B shows separating sets of extracted features 310 associated with 
the complexes 110, 120 ... on a continuous basis from the continuously extracted 
features 320 shown in Figure 3A. In the example embodiment shown in Figure 
3B, a set of features 340 associated with first complex is separated upon 
detecting the second subsequent complex 120. Also shown are separating a set of 
features 350 associated with second complex upon detecting a third subsequent 
complex. This process of separating a set of features associated with a complex 
is done on a continuous basis. Features associated with first complex are 
separated from the features of subsequent second complex in order to (so that the 
features are not mingled with the extracted features of the second subsequent 
complex) classify the sensed first complex 110. The next step is to identify 
fiducial features 360, 370, and so on. In the example embodiment shown in 
Figure 3B, fiducial feature 360 is identified based on comparing times when the 
separated sets of extracted features associated with the first complex 110 occur 
with a time when the heart beat associated with the first complex 1 10 is detected 
('*/ ), and choosing a feature from the set of features 340 that is closest in time 
to the time when the heart beat associated with the first complex 1 10 is detected. 
In the example shown in Figure 1, the first complex 1 10 is detected at 1 1 
milliseconds. The times when the separated set of features 340 associated with 
the first complex 110 occur are 2, 4, 8, 10, 12, 14, 20, and 24 milliseconds as 
shown in Figure 3 A. Based on comparing the times when the separated set of 
features 340 occur with the time when first complex 1 10 was detected, the 
feature 360 that occurs at a time of 10 milliseconds is the closest in time to the 
detected first complex 110 which occurs at 11 milliseconds. Therefore, extracted 
feature 360 is selected as the fiducial feature from the separated set of features 
340 associated with the first complex 1 10. Similarly, the extracted feature 370 is 
selected as the fiducial feature for the separated set of features 350 associated 
with subsequent second complex 120. This process repeats itself continuously on 
a real time basis. In some embodiments, the fiducial feature is identified from the 
set of separated features associated with the first complex based on a 
predetermined deviation value. The predetermined deviation value is based on a 
sample point having an amplitude farthest from a predetermined reference point. 
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The next step is to align the separated set of features 340 associated with 
the first complex 1 10, and to classify the separated set of features 340 by 
comparing the aligned set of features with a set of templates. Separated sets of 
features 340 and 350 are aligned to normalize the separated sets of feature 340 
and 350 so that the aligned features aid in comparing the separated sets of 
features with a template to classify the sensed complexes 1 10 and 120. Figure 
3B illustrates one example embodiment of aligning the set of extracted features 
340 (associated with first complex 1 10) with respect to an identified fiducial 
feature 360. In this example embodiment, the time when the identified fiducial 
feature occurs (at 10 milliseconds) is subtracted from each of the times when the 
separated set of features 340 associated with the first complex 110 occurs (at 2, 
4, 8, 10, ... milliseconds) to arrive at aligned times -8, -6, -2, 0, 2, 4, 10, and 14 
milliseconds for the separated set of features 340. Similarly, the aligned times -8, 
-2, 0, 2, 4, and 10 are computed for the separated set of features 350 associated 
with the second complex 120. This process of aligning separated set of features 
continues on a real time basis for the detected complexes on the sensed cardiac 
signal 100. It can be seen from the above illustrated computation, that the 
process aligning helps normalize the times associated with each of the separated 
sets of features, and how it aids in comparing the separated set of features 340 
and 350 with a set of templates and in classifying the sensed complexes 110 and 
120. 

Figure 4 illustrates a method 400 of identifying and classifying a sensed 
cardiac signal on a real-time basis. Method 400, as shown in Figure 4, begins 
with step 410 by sensing a cardiac signal 410 from the one or more electrodes 
disposed in or around a heart 515 on a continuous basis. 

The next step 420 in the process comprises computing curvatures at 

sample points X u X 2 , X 3> X x on the sensed cardiac signal. In some 

embodiments, the sample points X l5 X^ X 3> X x are fixed intervals of time at 

which the sensed cardiac signal is computed for its curvature. The sensed cardiac 
signal can be an atrial or ventricular signal of a heart. In one embodiment, the 

curvatures at sample point X la X» X 3 ^ X 1 is computed by fitting a cubic 

least square error curve to the sensed cardiac signal by using N sample points. In 
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this embodiment, N is an odd number greater than or equal to 5. Also in this 
embodiment, the cardiac rhythm management system computes the curvature at 
a mid point of the N sample points. In one embodiment, the cardiac rhythm 
management system computes curvature (K) at a sample point X x on the cardiac 
signal by using 5 sample points over which the cubic fit is (N being 5 in this 
example) made using 

K = (2C I /(1+B I 2 ) 3/2 ) 

where 

B, = T X Y(L2) + T 2 Y(M) + T 3 Y(I) +.T 4 Y(I+1) + T 5 Y(I+2) 

C, = S,Y(I-2) + S 2 Y(I-1) + S 3 Y(I) + S 4 Y(M) + S 5 Y(I+2) 

where S„ S 2 , S 3 ... and T„ T 2 , T 3 ...are constant terms in the summation. In this 
embodiment, I corresponds to the sample point at which the curvature is being 
computed. For example, if the sample point is X 3 , then I will be equal to 3 and so 
on. 

In another embodiment, the system computes average curvature between 
two adjacent sample points by using a linear interpolation of the B t and Q 
between the two adjacent sample points (this is done by doing the integral of the 
curvature between the two adjacent sample points). This integration helps 
capture any curvatures that become larger between the two adjacent sample 
points, so that the point-wise curvature values may more robustly represent the 
curvatures in the sensed cardiac signal. The system also integrates curvature 
versus time curve to find the area under the curve. This area under the curve 
represent the total change in angle between two points on the original cardiac 
signal. Thus, peaks, turns, inflections, and the like in the cardiac signals can be 
identified when the curvature versus time curve has peaks or nadirs. These peaks 
and nadirs denote points along the sensed cardiac signal where the turns are 
locally sharpest. Also, slow but steady turns on the sensed cardiac signal can be 
characterized by a small peak having a substantial area under the curvature 
curve. Thus, integrating curvatures between two adjacent sample points can help 
capture any curvatures that may become larger between any two adjacent sample 
points. 
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In one embodiment, the cardiac rhythm management system computes 

curvatures at sample points X„ X 2 , X 3 \ u X } on the sensed cardiac signal on a 

continuous basis. 

Then the next step 430 in the process includes extracting features from 
the computed curvatures by comparing the computed curvatures with a set of 
predetermined threshold values. Figure 3 A, illustrates extracting features from 

the computed curvatures at the sample points X„ X 2 , X 3 _ X L In some 

embodiments, the set of predetermined threshold values are based on a previous 
curvature value, a first curvature value, and a curvature threshold limit. In one 
embodiment, the step of extracting features from the curvature versus time signal 
is done on a continuous basis on the sensed cardiac signal as shown in Figure 
3 A. In some embodiments, the features are defined by one or more metrics. In 
this embodiment, the one or more metrics are defined as area under a computed 
curvature, a time of centroid of the area, and a value of original signal amplitude 
at a time of centroid of the area associated with the feature. 

Then the next step 440 comprises separating a set of features associated 
with a first complex upon detecting a second subsequent complex from the 
sensed cardiac signal. Separating of extracted features associated with the 
detected complex is illustrated in detail Figures 1, 2, 3 A and 3B. In some 
embodiments, the second subsequent complex is detected based on continuously 
processing the sensed cardiac signal to produce a signal that is the absolute value 
of the first derivative of a sensed cardiac signal versus time. Then comparing the 
produced absolute value versus time signal with a predetermined decaying 
threshold value to detect an occurrence of the second subsequent complex. In 
some other embodiments, the second subsequent complex is detected when the 
absolute value of the first derivative exceeds a predetermined decaying threshold 
value. In some embodiments, the features associated with the first complex are 
stored in a memory until the second subsequent complex is detected, and at that 
point, the features associated with the first complex are replaced by the features 
associated with the second subsequent complex in the memory. The step of 
separating can comprise separating the set of features based on selecting features 
associated with the first complex having a predetermined time earlier than the 
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detected second subsequent complex. In some embodiments, each of the set of 
features are further compared to an area threshold value to eliminate features less 
than the area threshold value. This can result in a set of revised features that can 
represent the first complex. 

Then the next step 450 in the process is to identify a fiducial feature from 
the set of separated features associated with the first complex based on a 
predetermined deviation value. In some embodiments, fiducial feature is 
identified based on comparing the time when the first complex occurs with the 
times when the separated set of features associated with the first complex occur. 
The process of identifying a fiducial feature is discussed in more detail with 
respect to Figures 3 A and 3B. 

Then the next step 460 in the process comprises aligning the set of 
separated features with respect to the identified fiducial feature. The process of 
aligning is also discussed in more detail in Figure 3B. The process of aligning 
helps normalize the times associated with each of the separated set of features 
and further facilitates in comparing the separated set of features with a set of 
predetermined templates to classify the detected complex. 

Then the next step 470 comprises comparing the aligned set of separated 
features with a set of predetermined templates. In one embodiment, the set of 
predetermined templates are one or more known types of heart beat signals. In 
another embodiment, the set of predetermined templates consists of sets of 
template boxes or zones for each of the known types of heart beat signals. 
Generally, each template box or zone can have center amplitude, a time width, 
and an amplitude width. Also generally, the template boxes or zones for each 
known type of heart beat can consist of one or more template boxes or zones 
each of which can have different centers and widths. In this embodiment, the 
comparing step comprises assigning a score based on how well each feature 
matches with each of the set of template boxes or zones. In this embodiment, if 
the separated set of features associated with the first complex fails to match with 
the set of templates then the complex is assigned with an unknown heart beat 
classification. In one embodiment, the set of template boxes or zones are 
modified over time to reflect any changes in sets of separated features used in 
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matching with the template boxes or zones. This can permit known beat type 
templates to follow slow trends or long-term changes in the cardiac signal. 

The next step 480 in the process comprises classifying the detected first 
complex based on the outcome of the comparison. Then the next step 490 
comprises providing a therapy based on the outcome of the classification. In one 
embodiment, the next step includes guiding a therapy to a heart based on the 
outcome of the classification. In another embodiment, the next step in the 
process includes storing classifications for future diagnostic purposes. 

Figure 5 is a schematic/block diagram 500 illustrating one embodiment 
of portions of a cardiac rhythm management system 505, which is coupled to a 
heart 515. Cardiac rhythm management system 505 includes a power source 580, 
a controller 540, a sensing circuit 520, a therapy circuit 590, and a via node/bus 
530. The sensing circuit 520 is coupled by a lead 510 to the heart 515 for 
receiving, sensing, and or detecting electrical heart signals. The sensing circuit 
520 provides one or more sensed cardiac signals to the controller 540, via the 
node/bus 530. The controller 540 also controls the delivery of a therapy provided 
by the therapy circuit 590 and/or other circuits, as discussed below. 

The controller 540 includes various modules, which are implemented 
either in hardware or as one or more sequences of steps carried out on a 
microprocessor or other controller. Such modules are illustrated separately for 
conceptual clarity; it is understood that the various modules of the controller 540 
need not be separately embodied, but may be combined and/or otherwise 
implemented, such as in software/firmware. 

In general terms, the sensing circuit 520 senses a cardiac signal on a 
continuous basis from a heart tissue in contact with a catheter lead 5 10 to which 
the sensing circuit 520 is coupled. In one embodiment, the sensed cardiac signal 
is an atrial signal. In another embodiment, the sensed cardiac signal is a 
ventricular signal. Sensed cardiac signal from the sensing circuit 520 is then 
received and processed by an analyzer 550 of the controller 540 to compute 

curvatures continuously on a real-time basis at sample point X u X 2 , X 3 Xj 

on the sensed cardiac signal. In some embodiments, the sample points are spaced 
at fixed intervals of time preset by the analyzer 550 on the sensed cardiac signal. 
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In one embodiment, the analyzer 550 computes curvature at sample points X„ 
X 2 , X 3 X, by fitting a cubic least square error curve to the sensed cardiac 
signal using N sample points to fit the cubic least square error curve. In one 
embodiment, the analyzer 550 computes the computing curvature (K) at a 
sample point X, on the cardiac signal when using 5 sample points over which the 
cubic fit is (N being 5 in this example) made using 
K=(2C,/(1+B, 2 ) 3/2 ) 

where 

B, = T,Y(I-2) + T 2 Y(I-1) + T 3 Y(I) + T 4 Y(I+1) + T s Y(I+2) 

C, = S,Y(I-2) + S 2 Y(I-1) + S 3 Y(I) + S 4 Y(M) + S s Y(I+2) 

where S„ S 2> S 3 ... and T„ T 2 , T 3 ...are constant terms in the summation. In this 
example embodiment, I is equal to the point at which the curvature is being 
computed. For example, if the curvature is being computed at sample point X 3 , 
then I will be equal to 3 and so on. In one embodiment, N is an odd number 
greater than or equal to 5. In this embodiment, the analyzer 550 computes 
curvature at a midpoint of the N number of sample points. 

In another embodiment, the analyzer 550 computes curvature based on an 
average curvature between two adjacent sample points based on a linear 
interpolation of the B, and Cj between two adjacent points (integrating the 
computed curvatures between the two adjacent sample points). In one 
embodiment, the analyzer 550 computes the curvatures continuously on a real- 
time basis from the sensed cardiac signal. Then the analyzer 550 extracts features 
by comparing the computed curvatures to a set of predetermined threshold 
values. In one embodiment, the predetermined threshold values are based on a 
previous curvature value, a curvature value, and a curvature threshold limit. In 
one embodiment, the analyzer 550 extracts features on a continuous basis. Then 
the analyzer 550 separates a set of features associated with a first complex upon 
detecting a second subsequent complex from the sensed cardiac signal. Then the 
analyzer 550 identifies a fiducial feature from the set of separated features 
associated with the first complex based on a predetermined deviation value. 
Then the analyzer 550 further aligns the set of separated features with respect to 
the identified fiducial feature. In one embodiment, the analyzer 550 separates the 
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set of features associated with the first complex based on a predetermined time 
earlier than the detected second subsequent complex. In one embodiment, the 
predetermined deviation value is based on a sample point having an amplitude 
farthest from a predetermined reference point. 

Then a comparator 570 coupled to the analyzer 550, compares the 
aligned set of features associated with the first complex with a predetermined set 
of templates, and classifies the first complex based on the outcome of the 
comparison. In one embodiment, the comparator 570 issues a command signal 
based on the outcome of the classification. In another embodiment, the therapy 
circuit 590 coupled to the comparator 570 delivers the electrical energy through 
the lead 510 to at least one electrode disposed in or around the heart 515 upon 
receiving the command signal from the comparator 570. The electrical energy 
delivered by the therapy circuit 590 can be a pacing pulse electrical energy, hi 
one embodiment, the analyzer 550 includes a variable gain circuitry that can 
adapt according to the changes in the sensed cardiac signal. In one embodiment, 
the controller 540 includes a memory 560 coupled to the comparator 570 and the 
analyzer 550 to store the extracted features of the first complex. In another 
embodiment, the memory 560 stores the classified first complexes for diagnostic 
purposes. 

It is to be understood that the above description is intended to be 
illustrative, and not restrictive. Many other embodiments will be apparent to 
those of skill in the art upon reviewing the above description. The scope of the 
invention should, therefore, be determined with reference to the appended 
claims, along with the full scope of equivalents to which such claims are entitled. 

Conclusion 

The above described method provides, among other things, a curvature 
based complex identification and classification. The process comprises sensing a 
cardiac signal and computing curvatures at sample points on the sensed cardiac 
signal, extracting features from the computed curvatures, comparing the 
extracted features with a set of predetermined templates, and classifying the 
cardiac signal based on the outcome of the comparison. 



15 



WO 02/40094 



PCT/US01/45763 



WHAT TS CLAIMED IS: 

1 . A cardiac rhythm management system, comprising: 

at least one electrode; 
5 a signal sensing circuit coupled to the electrode to sense a cardiac signal; 

a controller coupled to the sensing circuit, wherein the controller receives 
the sensed cardiac signal, and wherein the controller includes: 

an analyzer, to compute curvatures at sample points X t , X 2 , X 3 X! 

on the sensed cardiac signal, wherein the analyzer extracts features by comparing 
10 the computed curvatures to a set of predetermined threshold values; and 

a comparator, coupled to the analyzer, compares the extracted features 
with a set of predetermined templates, and classifies the sensed cardiac signal 
based on the outcome of the comparison. 

15 2. The system of claim 1 , wherein the sensing the cardiac signal includes 
sensing complexes on a real-time basis, wherein the complexes comprise heart 
beat signals. 

3 . The system of any of the preceding claims, wherein computing 
20 curvatures comprises computing curvatures at the sample points X„ X 2 , X 3f 

Xj by fitting a cubic least square error curve to a first complex of the sensed 
cardiac signal using an N sample points to fit the cubic least square error curve, 
wherein the analyzer extracts features of the first complex by comparing the 
computed curvatures of the first complex to a set of predetermined threshold 

25 values, and wherein the analyzer further separates a set of features associated 
with the first complex upon detecting a second subsequent complex from the 
sensed cardiac signal; and wherein the comparator compares the aligned set of 
features associated with the second subsequent complex with a set of 
predetermined templates, and classifies the first complex based on the outcome 

30 of the comparison. 
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4. The system of any of the preceding claims, wherein the comparator 
issues a command signal based on the outcome of the classification. 

5. The system of any of the preceding claims further comprises a therapy 
5 circuit, coupled to the comparator, to deliver electrical energy through the at 

least one electrode upon receiving the command signal from the comparator. 

6. The system of claim 5, wherein the electrical energy is a pacing pulse 
electrical energy. 

10 

7. The system of claim 3, wherein the set of predetermined threshold values 
are based on a previous curvature value, a curvature value, and a curvature 
threshold limit 

15 8. The system of claim 3, wherein the analyzer further identifies a fiducial 
feature from the set of separated features associated with the first complex based 
on a predetermined deviation value, and further aligns the set of separated 
features associated with the first complex with.respect to the identified fiducial 
feature. 

20 

9. The system of claim 8, wherein the analyzer separates the set of features 
associated with a first complex based on a predetermined time earlier than the 
detected second subsequent complex. 

25 10. The system of claim 9, wherein the predetermined deviation value is 
based on a sample point having an amplitude farthest from a predetermined 
reference point. 

11. The system of claim 3, wherein N is an odd number greater than or equal 
30 to 5. 
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12. The system of claim 1 1, wherein the computing curvature comprises 
computing curvature at a mid point of the N number of sample points. 

13. The system of claim 12, wherein the analyzer computes the curvature (K) 
5 at a sample point X x on the cardiac signal when using 5 sample points to fit the 

cubic least square error curve, based on 



K = (2C I /(1+B 1 2 ) 3/2 ) 

where 

10 B, = ^¥(1-2) + T 2 Y(I-D + T 3 Y(I) +.T 4 Y(I+ 1) + T 5 Y(I+2) 
Q - + S 2 Y(I-1) + S 3 Y(I) + S 4 Y(I+1) + S 5 Y(I+2) 

where S and T are constants. 



14. The system of claim 1 3, wherein the analyzer computes an average 

1 5 curvature between two adjacent sample points based on linear interpolation of 
theB! and Q between two adjacent points and integrating the computed 
curvatures between the two adjacent sample points. 

15. The system of any of the preceding claims, wherein the analyzer further 
20 comprises a variable gain, wherein the variable gain adapts according to changes 

in the sensed cardiac signal. 

16. The system of any of the preceding claims, wherein the at least one 
electrode is disposed in or around a heart. 

25 

1 7. The system of any of the preceding claims further comprises a memory to 
store the extracted features of the first complex. 

1 8. The system of claim 1 7, wherein the memory further stores the classified 
30 first complexes for diagnostic purposes. 

19. A method, comprising: 
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sensing a cardiac signal; 

computing curvatures at sample points X u X 2 , X 3 X l on the sensed 

cardiac signal; 

extracting features from the computed curvatures; 
5 comparing the extracted features with a set of predetermined templates; 

and 

classifying the cardiac signal based on an outcome of the comparison. 

20. The method of claim 19, wherein the cardiac signal comprises sensing 
1 0 complexes in real time, wherein the complexes are cardiac cycles. 

2 1 . The method of any of the preceding claims, wherein computing the 

curvatures at sample points Xj, X 2 , X 3 X { comprises computing curvatures 

at the sample points X lf X 2 , X 3j X r by fitting a cubic least square error 

1 5 curve to sensed complexes on the cardiac signal, wherein the sample points have 
a predetermined interval of time between adjacent sample points. 

22. The method of claim 21, wherein computing the curvatures at sample 
points X„ X 2 , X 3> X r comprises fitting a cubic least square error curve 

20 around a sample point by using N number of sample points to fit the cubic least 
square error curve, wherein N is an odd number greater than or equal to. 5. 

23. The method of claim 22, wherein computing the curvature at the sample 
point using the N number of sample points comprises using the sample point as a 

25 mid point of the N number of sample points. 

24. The method of claim 23, wherein computing the curvature (K) at a 
sample point X, of the cardiac signal when using 5 sample points to fit the cubic 
least square error curve, is based on 

30 

K = (2C ! /(1+B I 2 ) 3 ' 2 ) 

where 
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B, = T,Y(I-2) + T 2 Y(I-1) + T 3 Y(0 +.T 4 Y(M) + T 5 Y(I+2) 

C, = S^Ct-2) + SjYCI-l) + S 3 Y© + S 4 Y(I+1) + S 5 Y(I+2) 
where S and T are constants. 

5 25. The method of claim 24, wherein computing the curvature comprises 
computing an average curvature between two adjacent sample points based on a 
linear interpolation of the Bj and C, between two adjacent points and integrating 
the computed curvatures between the two adjacent sample points X x> X 2 , 

X P 

10 

26. The method of any of the preceding claims, wherein extracting features 
from the computed curvatures includes: 

comparing the computed curvatures to a set of predetermined threshold 
values; and 

1 5 identifying and separating a set of extracted features associated with a 

first complex upon detecting a second subsequent complex from the sensed 
cardiac signal. 

27. The method of claim 26, wherein the set of predetermined threshold 
20 values are based on a previous curvature value, a first curvature value, and a 

curvature threshold limit. 

28 . The method of claim 27, wherein separating the set of extracted features 
comprises separating the set of features based on identifying features associated 

25 with the first complex having a predetermined time earlier than the detected 
second subsequent complex. 

29. The method of any of the preceding claims, wherein comparing the 
extracted features comprises comparing the separated set of extracted features 

30 associated with the first complex with a set of predetermined templates. 
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30. The method of any of the preceding claims, wherein comparing the 
extracted features further comprises: 

identifying a fiducial feature from the set of separated features associated 
with the first complex based on a predetermined deviation value; and 
5 aligning the set of separated features associated using the identified 

fiducial feature. 

3 1 . The method of claim 30, wherein the predetermined deviation value is 
based on a sample point having an amplitude farthest from a predetermined 

10 reference point. 

32. The method of any of the preceding claims further includes repeating the 
above steps for a real-time classification of heart beat signals from the sensed 
cardiac signal. 

15 

33. The method of claim 31, wherein the feature is defined by one or more 
metrics. 

34. The method of claim 33, wherein the one or more metrics are area under 
20 a computed curvature, a time of centroid of the area, and a value of original 

signal amplitude at a time of the centroid of the area. 

35. The method of claim 31, wherein comparing the set of features comprises 
comparing the set of features associated with the first complex with one or more 

25 predetermined heart beat signals. 

36. The method of claim 26, wherein the predetermined set of templates 
comprises one or more predetermined template zones defined by a center time, a 
center amplitude, a time width, and an amplitude width. 

30 

37. The method of any of the preceding claims further comprises providing a 
therapy to a heart based on the outcome of the classification. 
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38. The method of any of the preceding claims further comprises guiding a 
therapy to a heart based on the outcome of the classification.' 

39. The method of any of the preceding claims further comprises storing 
5 classifications for diagnostic purposes. 

40. The method of any of the preceding claims, wherein computing the 
curvature further comprises permitting the computed curvature signal to have a 
variable gain that adapts according to the changes in sensed cardiac signal. 
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